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ABSTRACT 

Context. The intensity of the Ca II K resonance line observed with spectrographs and Lyot-type filters has long served 
as a diagnostic of the solar chromosphere. However, the literature contains a relative lack of photometric measurements 
of solar features observed at this spectral range. 

Aims. We study the radiative emission of various types of solar features, such as quiet Sun, enhanced network, plage, 
and bright plage regions, identified on filtergrams taken in the Ca II K line. 

Methods. We analysed full-disk images obtained with the PSPT, by using three interference filters that sample the 
Ca II K line with different bandpasses. We studied the dependence of the radiative emission of disk features on the filter 
bandpass. We also performed a non-local thermal equilibrium (NLTE) spectral synthesis of the Ca II K line integrated 
over the bandpass of PSPT filters. The synthesis was carried out by utilizing the partial frequency redistribution (PRD) 
with the most recent set of semi-empirical atmosphere models in the literature and some earlier atmosphere models. 
As the studied models were computed by assuming the complete redistribution formalism (CRD), we also performed 
simulations with this approximation for comparison. 

Results. We measured the center-to-limb variation of intensity values for various solar features identified on PSPT 
images and compared the results obtained with those derived from the synthesis. We find that CRD calculations derived 
using the most recent quiet Sun model, on average, reproduce the measured values of the quiet Sun regions slightly 
more accurately than PRD computations with the same model. This may reflect that the utilized atmospheric model 
was computed assuming CRD. Calculations with PRD on earlier quiet Sun model atmospheres reproduce measured 
quantities with a similar accuracy as to that achieved here by applying CRD to the recent model. We also find that the 
median contrast values measured for most of the identified bright features, disk positions, and filter bandpasses are, 
on average, a factor ~ 1.9 lower than those derived from PRD simulations performed using the recent bright feature 
models. The discrepancy between measured and modeled values decreases by « 12% after taking into account straylight 
effects on PSPT images. When moving towards the limb, PRD computations display closer agreement with the data 
than performed in CRD. Moreover, PRD computations on either the most recent or the earlier atmosphere models of 
bright features reproduce measurements from plage and bright plage regions with a similar accuracy. 

Key words. Sun: activity - Sun: photosphere - Sun: chromosphere - Sun: faculae, plages 



. 1. Introduction us to check the various semi-emp irical models used 
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^ ■ Observations taken with spectrographs and Lyot-type fil- | Reardon . Uitenbroek fc CauzzillM iPietarila et al.l l20091 

H , ters m the Ca II K resonance line have l ong served as and the recon structi on of irradia nce v ariations (e.g. 

^ , diagnostics of the solar chromosphere (e.g., | Rutterj| | 2007t | Domhlg0j et al.l [20091 lHarder et all l200l iKrivova et all 

and references therein) and have been used to deter- | 2009h In part icular. mo dels have been presented 

^^^ f^J^^ ^ ^S^ ^ Vontenla et al-iH that differ from previous sets of mod- 



sphere mode s fenazza , etoajj [1981J | Maltby et ajj [19861 ds m terms of the physical assumpt ions used to derive their 
IFontenla et all 119911 119931). However, exc ept for the quiet tpmnprahlrp structure These models have been adapted in 
Sun ^rkei|1968; White Si Wittmannfel 
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is a shortage of measurements m the literature of the ra- } th ^ com ison with observations in the near UV 
diative emission of solar features observed m the Ca II K range wag Umited Detailed comparisons betwe e n measure- 
range at various heliocentric angles. ments and model predictions in the near rjy^ inc i uding for 
In addition to providing new data, photometric the c& n K lin6j WQuld aid in model refinement . 

measurements of disk features in Ca II K allow . ,,, . r . .. ,. 
In this paper, we present a study of the radiative emis- 
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Fig. 1. Transmission profiles of the three PSPT interfer- 
ence filters considered in this study. The three profiles, each 
normalized to its maximum, are represented with different 
colors and lines as indicated in the legend. They are super- 
imposed on the solar reference spectrum (black line) pro- 
vided bv lDelbouille et all (jl972ft . 



bright plage regions, identified on full-disk observations 
taken in the spectral range of the Ca II K line. The images 
were obtained with the PSPT (Precision Solar Photometric 
Telescope), by using three interference filters that sample 
the Ca II K line with different bandpasses. We studied the 
dependence of the measured emission on the filter band- 
pass. We also compared the results obtained with those de- 
rived from the numerical synthesis perfor med for the band- 
pass of PSPT filters with the RH c ode flUitenbroeld |2002) 
utilizing the new models presented bv lFontenla et all ([2 009) 
as well as earlier atmosphe re models ([Vernazza et al.|[l981 ; 
iFontenla et "al~lll993l 12006ft . 



2. Observations and data reduction 

The PSPT is a 15 cm, low-scattered-light, refracting tele- 
scope designed for synoptic photometric solar observa- 
tions characterized by 0.1% pixel-to- p ixel relative pho- 
tometric precision (|Coulter et all Il994h . Two such tele- 
scopes, which differ only in minor hardware characteris- 
tics and operational strategies, are operated at Mauna 
Loa (MLSO-PSPT) and Monte Porzio Catone (Rome- 
PSPT) Observat ories by the High Altitude Observatory 
Ras t et alj 11999ft an d Osservatorio Astronomico di Roma 



Ermolli et al.lll998ft . respectively. These telescopes typi- 



cally acquire full-disk solar images on 2048x2048 CCD ar- 
rays with narrow-band interference filters centred on the 
blue continuum (409.412 FWHM 0.267 nm), red continuum 
(607.095 FWHM 0.458 nm), and Ca II K (393.415 FWHM 
0.273 nm, hereafter PK-027). However, the installation of 
additional filters on both the telescopes has also allowed the 
acquisition of images at other spectral ranges. In particu- 
lar, two interference filters sampling the Ca II K range were 
added to the MLSO telescope in summer 2007, one at line 
center (393.379 FWHM 0.103 nm, hereafter PK-010) and 
the other in the red wing of the line (393.633 FWHM 0.106 
nm, hereafter PKR-011). These two filters were operated 
simultaneously with the broader band PK-027 filter. Here 
we analyse the observations carried out with these three 
Ca II K filters, hereafter referred to as the PK filters. 



Figure [T] shows the transmission profiles measured for 
the PK filters superimposed on a reference solar spectrum. 
The central wavelength and FWHM of these profiles are 
summarized in Table [21 The filter band widths of PK do not 
allow the features of the Ca II K line center to be resolved. 
These features, the core, the reversal (emission peaks), and 
the secondary minima occur within a spectral range less 
than 0.1 nm wide. 

We analysed 157 sets of observations available from the 
MLSO PSPT archivcQ for the period 2007 June 7 to July 
31. Each set is composed of the images obtained with the 
PK filters and one in the red continuum. The red continuum 
images are, on average, separated in time from the Phones 
by 7-8 minutes. 

We separated the data into two samples, corresponding 
to the highest and the lowest quality images. They con- 
tain 37 and 33 sets of images, respectively. Image qual- 
ity was determined by measuring the width of the solar 
limb in the PK-027 images, which provides an estimate 
of the atmospheric seeing cond itions (jCriscuoli et alJl2007t 
iRast. Ortiz and Meisnerll2008ft . In particular, images with 
an average limb width smaller than 4 pixels constitute the 
highest quality (hereafter referred to as good) images, while 
the ones with average limb width larger than 5.5 pixels form 
the set of the worst (hereafter referred to as poor) images. 
The average and the standard deviation a of limb width 
values measured for the whole sample of PK-027 images 
are 4.7 ± 0.8 pixels; the thresholds introduced above are 
the values derived from the average ±<j. 

All the images analysed for our study were pre- 
processed to apply the correction for gain variations across 
the CCD device and to compensate for residual large-scale 
linear gradients affecting the im ages after CCD calibration 
(|Rast. Ortiz an d Meisner 2008). The images were re-sized 
and aligned, by using linear interpolation, to a common ref- 
erence grid on which the solar disk size and orientation are 
constant. 

During the period analysed, the Sun was in a quiet 
phase. However, eight small active regions appeared on the 
solar disk, specifically AR 10958 to AR 10965, in addition to 
some plages and enhanced network regions. To separate the 
various disk features seen in PK ima ges, we used th e imag e 
decomposition method described by Fontenla et all ([2009). 
This method utilizes a threshold scheme derived from the 
partitioning of intensity histograms constructed from the 
images as a function of heliocentric angle (/z) . The scheme 
assigns structure types based on the agreement between the 
intensity measured at each image pixel and the one defined 
by the histogram partition. In other words, features are de- 
fined by their normalized intensity as a function of /z, and 
pixels arc identified with particular features based on their 
observed position and intensity. 

The decomposition code used in our study assumes 
seven classes of disk features identified in PK-027 and red 
continuum co-aligned images. Following b oth the scheme 
and r eference atmospheres proposed by IFontenla et al.l 
()2009ft . these classes are coded by the letters B, D, F, H, 
P, S, and R. Features B through P represent the intensity 
criterion applied to pixels of PK-027 images. They local- 



1 Calibrated images from the synoptic observations 
taken with the tw o PSPTs over many years are avail- 
able on-line at http://www.oa-roma.inaf.it/solare and 
http://lasp.colorado.edu/pspt_ access/ . 
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Table 1. Identification classes, solar features, and atmo- 
sphere models utilized in this study. 



Class 


Solar Feature 


Atmosphere model 


B 


quiet Sun 


model B, VAL3-C, FAL3-C, FA-06 


D 


network 


model D, FAL3-F 


F 


enhanced 


model F 




network 




H 


plage 


model H, model F, FAL3-F, FAL3-P 


P 


bright plage 


model P 


S 


umbra 




R 


penumbra 





Notes. The table summarizes the classes utilized for the iden- 
tification of the various solar features, the relation between each 
class and the corresponding solar feature, as well as the relation 
between solar features and atmosphere models considered for the 
comparisons presented in the following. Models B to P indicate 
the set of atmospheres of Fontenla et al. (2009), which consti- 
tutes the reference set in this study. Models VAL3-C, FAL3- 
C, and FA06-C indicate the atmospheres for quiet Sun regions 
presented by Vernazza et al. (1981), and Fontenla et al. (1993, 
2006), respectively. Models FAL3-F and FAL3-P are the atmo- 
spheres for plage and bright plage regions presented by Fontenla 
et al. (1993). 



ize average median quiet Sun (internetwork, QS hereafter), 
network (QS network lanes), enhanced network, plage, and 
bright plage (or facula) regions, respectively. Features S 
and R, which correspond to umbral and penumbral regions, 
respectively, are identified in red continuum images. The 
threshold values utilized in PK-027 images are shown in 
the top panels of Figs. [5] and [7] We note that feature B 
is called here quiet Sun (QS), but in the threshold scheme 
localizes more broadly internetwork regions, which likely 
contain residual magnetism. We also note that the mean- 
ing of quiet Sun here is quite distinct from that often used 
in describing chromospheric or transition region emission, 
where it is often used to mean the non-coronal hole por- 
tion of the Sun. Table [T] summarizes the classes utilized 
to identify the various solar features, the relation between 
each class and the corresponding solar feature, as well as 
the relation between solar features and the various atmo- 
sphere models considered in the comparisons presented in 
the following. 

Figure shows examples of features identified in two 
sets of observations analysed in this study. We measured the 
radiative emission of the features identified in each image. 
In particular, we considered the emergent intensity of quiet 
Sun regions and the contrast of bright features. The latter 
quantity is defined, for each heliocentric angle, as the ratio 
(If — Ib)/Ib, where // is the median of intensity values of 
pixels labeled as feature /, and Ib is the median intensity of 
pixels coded as feature B (QS, internetwork). The contrast 
of features was computed in 50 equal-area annuli centered 
on the solar disk center. For each bandpass and feature, 
the center-to-limb variation (CLV, hereafter) of contrast 
was computed by taking into account the median contrast 
values measured for the given feature in the 50 annuli of all 
the images in the specified bandpass. We note that we de- 
composed PK-027 images, and subsequently measured the 
radiative properties of identified features in PK-027, PKR- 
011, and PK-010 co- aligned images using this decomposi- 
tion. This was done to compare measurement and synthesis 
results assuming that the size of identified features is con- 




Fig. 2. Example of PK-027 (top panels) observations anal- 
ysed in this study and of the corresponding mask images 
(bottom panels). The observations were taken on 2007 June 
7 (left side) and June 21 (right side). The disk features on 
mask images are identified by color: blue, light blue, green, 
and yellow show quiet Sun, network, enhanced network and 
plage regions, respectively; white and red show penumbral 
and umbral regions (visible only when enlarging the figure) . 



stant across the various images, i.e., over the atmospheric 
heights sampled by the PK filters. 

3. Spectral synthesis 

To discuss the results of our measurements and place them 
in context with results of previous studies, we computed 
synthetic spectra in th e range of the C a II K line with 
the RH radiative code ([Uitcnbrock 120021) . This code per- 
forms NLTE radiative transfer modeling by utilizing the 
partial frequency redistribution formalism (PRD). It has 
been utilized for spectral synthesis and atmosphere di- 
agnostics at chromospheric heigh ts employing var i ous at - 
mosphere models. For examp l e. iTritschler et all (l2007f), 
iReardon. Uitenbroek fc Cauzzil (|2009[) . and iPietarila et aL 
(2009) utilized RH and ear l ier mo dels by iFontenla et aL 
jl993h . while iRezaei et all (120081) and iGrigorveva et al. 
(200 9j) employ e d the sets of IFontenla et al.l (|2006) and 
IFontenla etail (|2007f ). 

Our calculations with RH were perfo rmed by adopting 
a 5-le vel plus continuum atomic model (lUitenbroekl ll989. 
I2001D . the atomic data of IShine fc Linskvl (119741) wit h col- 
lisional strength values from iMelendez et al. | (120071). and 
the se t of atmosphere models presented bv IFontenla et al.l 
(120091) . hereaft er FC0 9, as well as the ear li er models o f 
IVernazza et"aU f|l98lh and IFontenla et"al1 (|l993l . l2006h . 
hereafter VAL3, FAL3, and FA06, respectively. Since some 
of these models were constructed with radiative transfer 
codes that utilize the complete redistribution approxima- 
tion (CRD), which is also an option in RH, we computed 
the Ca II K line profiles with both formalisms for compar- 
ison. This was done to evaluate the sensitivity of synthesis 
results to code approximations. 

For each reference model, we calculated the emergent 
spectrum at 43 wavelengths over the range 392.515 nm to 
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Fig. 3. Intensity response functions to perturbations of 
temperature for the PSPT filters computed at the disk po- 
sitions fi = 1 (top) and jtx = 0.4 (bottom). The functions 
derived for the three filters and various atmosphere models 
are indicated by different colors and lines as specified in the 
legend. Each function is normalized to its maximum. 



394.214 nm for 13 different viewing angles. For each an- 
gle and atmosphere, we weighted the obtained spectrum 
with the various filter profiles considered in our study and 
summed up the result over the spectral range of our synthe- 
sis. The intensities derived with the code were then com- 
pared with those measured in PK images for various disk 
features. The results obtained are presen ted in Sect. 4. 

As discussed by I Shine et"al~l ()1975| ) and lUitenbro'ekl 
l|2002h . we found that the Ca II K line profiles derived 
from our simulations using the various reference models 
and adopting either PRD or CRD, differ in terms of the line 
relative minima. The amplitude of the line emission peaks 
derived for the various models at disk center, however, is 
almost unaffected by the adopted redistri bution formal- 
ism. T his is in contrast to the suggestion by iFontenla et al.l 
(2009, their Fig. 7) that the much higher emission peaks 
in the H and K lines they found can be attributed to the 
lack of inclusion of PRD effects. We note that the differ- 
ence in emission peak amplitude betwee n CRD and PRD 
increases considerably towards the limb (|Uitenbroeki ri989) 
in any given model, and that this affects the CLV of line- 
integrated quantities. We also found that the difference be- 
tween the profiles derived from the two approximations in- 
creases according to the brightness of the model disk fea- 
ture it represents. Comparison of emission peaks obtained 
with both PRD or CRD also indicates that model P differs 
from the other models in terms of the smaller value of its 
micro-turbulent velocity. 

Next, we computed the intensity response function (RF) 
to perturbations of temperature for the PK. This func- 
tion (RFs-PK, hereafter) allows us to infer the atmospheric 
heights from the intensity m easured in PK ima ges more re- 
liably than other methods (|Uitenbroekl 120061 ) . Moreover, 
it also indicates the sensitivity of synthesized quantities 
to atmospheric models, though the CRD and PRD cal- 
culations are affected by uncertainties associated with the 
static and one-dimensional nature of the atmospheres used. 
We computed the RFs-PK only with PRD because of the 
advantages that t his approach offer s for the synthesis of 
the Ca II K line ([Uitenbroekl 120021 ) over the PK range, 



Table 2. Characteristics of the PK filters. 







PK-010 


PK-027 


PKR-011 


Ao 
BW 


(nm) 
(nm) 


393.379 
0.103 


393.415 
0.273 


393.633 
0.106 


/x— 1 
Ave 
FW 
RF-f 


(km) 
(km) 


375±25 
270±30 
0.58-0.82 


200±5 

260±15 

0.84-0.94 


175±5 
160±5 
0.996-0.998 


li=0A 
Ave 
FW 
RF-f 


(km) 
(km) 


380±30 
270±35 
0.43-0.76 


300±5 

260±15 

0.75-0.91 


200±5 
140±5 
0.989-0.995 



Notes. For each bandpass analysed, the table summarizes the 
central wavelength (Ao) and FWHM (BW) of the filter transmis- 
sion profile, the average response height (ave), and FWHM of 
response heights (FW) of the various response functions (RFs) 
derived from the atmosphere models of Fontenla et al. (2009) at 
disk positions fj,— l and fj,=QA, and the fraction of RF (RF-f) 
from atmospheric heights below 500 km at the same disk posi- 
tions. 



where the Ca II K dominates all the opacity and emissiv- 
ity sources. The filter-response functions of the different 
atmospheres were calculated numerically as described by 
lUitenbroekl (|2006i ). i.e. by perturbing the temperature of an 
atmosphere at each height, evaluating the emergent inten- 
sity integrated over the filter bandpass, and comparing the 
perturbed filter intensities with those of the un-perturbed 
atmosphere. 

Figure [3] shows the RFs-PK derived for the various at- 
mosphere models at the disk positions /i = 1 (top) and 
fi = 0.4 (bottom). These RFs-PK indicate that the PK sam- 
ple quite a wide range of atmospheric heights, and that the 
radiative signals measured with these filters are dominated 
by the wings of the Ca II K line, which form at heights 
below 500 km. We found that PK-010 samples the high- 
est atmospheric regions, while PKR-011 and PK-027 sam- 
ple middle and upper photospheric heights, respectively. In 
particular, PKR-011 has the sharpest RF and provides the 
cleanest photospheric signal. These properties follow from 
the filter widths and positions. 

All RFs-PK display a smooth variation with atmosphere 
height, except for the ones derived for model B (representa- 
tive of QS), which show a sharp variation at a height of « 
850 km produced by the temperature profile of this model 
above the temperature minimum layers. The RF derived for 
PK-010 and PK-027 also have a long tail extending to the 
middle and upper chromospheric heights exceeding 1000 
km. However, the response of PK-027 is k 3 times lower 
than that of PK-010 at both heights of 800 km and 1000 
km, for a vertical line of sight and the reference atmosphere 
model P (bright plage). In addition, while the RFs-PK de- 
rived for models B (QS) to model H (plage) exhibit similar 
variations with height, the RF-PK computed for model P 
(bright plage) is w 3 times higher than the RF-PK obtained 
for model H (plage) at a height of 800 km, but it is w 3 
times lower than the RF for model H (plage) at 1000 km 
height. 

Table [5] summarizes the average and FWHM of the re- 
sponse heights of the RFs-PK for the entire set of atmo- 
sphere models. We found that these quantities are only 
slightly sensitive to the choice of atmosphere models and fil- 
ter profiles. Among the various RFs-PK, the most sensitive 
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Fig. 4. Top: Quiet Sun limb darkening measured in PK-027 
(green, solid line), PK-010 (violet, dotted line), and PKR- 
01 1 (red, dashed line) good images. Intensity values are 
normalized to the median of values measured at the disk 
center (/i > 0.95). The error bars show the dispersion of 
measured values. Bottom: Ratio of the median of quiet Sun 
limb darkening measured in PK-010 to PK-027 (diamond) 
images, and in PKR-011 to PK-027 images (asterisk). The 
error bars show the dispersion of values from error propa- 
gation. The dashed line would represent perfect agreement 
between the two compared sets. 
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Fig. 5. Ratio between the median of quiet Sun limb dark- 
ening measured in PSPT images and the quiet Sun limb 
darkening derived from the synthesis performed on model 
B by FC09 for the PK-027 (green, solid line), PK-010 (vi- 
olet, dotted line), and PKR-011 (red, dashed line) niters 
by assuming either the PRD (top panel) or CRD (bottom 
panel) formalism. 



to chromospheric conditions are the ones derived for model 
P and PK-010. We also found that the average response 
height increases from 1% to 40% when moving toward the 
limb for PK-027 and PK-010. 



4. Results 

Figure H] (top panel) shows the CLV of intensity values 
measured for QS identified in PK images. In particular, 
the CLV computed from the PK-010 images is less steep 
than the one found in both PKR-011 and PK-027 obser- 



vations. Moreover, on average, more poorly resolved im- 
ages from the poor quality set lead to the CLV being 
larger than the median CLV measured for good quality im- 
ages, although differences are within the deviation of val- 
ues measured on good images. These findings are in agree- 
me nt with results obtained for d iffere nt Ca II K data sets 
bv iLivingston fc Sheelevl (|2008h and ICriscuoli k Ermollil 
(120081) . 

The dependence of the intensity CLV on both filter 
bandwidth and central wavelength is, however, small. We 
found that the average and the standard deviation of the 
ratio of the values measured on PK-010 to PK-027 images 
are 1.013 ± 0.013, for all the disk positions with fi > 0.2 
(Fig.Sl bottom panel). The same quantities derived for the 
ratio of the values measured on PKR-011 to PK- 027 images 
are 0.99 ± 0.01. We note that these averages are strongly 
weighted by the values obtained at positions close to disk 
center, because of the larger number of pixels corresponding 
to large \i. 

Comparing the median CLV measured in PK images 
with the CLV derived from line synthesis, we find that the 
measured CLVs are generally consistent with the results 
of the spectral synthesis performed on the reference atmo- 
sphere (model B) utilized in this study. In particular, Fig. [5] 
(top panel) shows this comparison for results derived with 
PRD. The differences between measured and modeled val- 
ues increase towards the limb, reaching between 4% and 
8% depending on the bandpass, with the measured limb 
values always brighter than model predictions. The average 
and the standard deviation of the ratio of the measured to 
modeled values for all the disk positions with >0.2 are 
0.999 ± 0.003, 1.015 ±0.015, and 1.013±0.016 for the PK- 
027, PK-010, and PKR-011 images, respectively. The same 
quantities evaluated for poor quality images are within the 
deviation of values measured in high quality data. These re- 
sults indicate that the intensity thresholds utilized in this 
study for the identification of QS are in fairly good agree- 
ment with those obtained with PRD in the reference QS 
atmosphere (model B) of FC09. Figure [5] (bottom panel) 
also shows that the agreement between measurements and 
synthesis results slightly improves by taking into account 
the results from CRD computations, in line with the in- 
crease in emission of line emission peaks towards the limb 
for profiles calculated with this formalism. We found that 
the maximum difference between measured and modeled 
values is ~ 5%. The average and standard deviation of the 
ratio of measured to modeled values obtained with CRD 
over all disk positions with fi >0.2 are 0.99 ± 0.01, 1.006 
±0.005, and 1.01±0.01 for the PK-027, PK-010, and PKR- 
011 images, respectively. We discuss these results in Sect. 
5. 

We now consider the contrast measured for the vari- 
ous bright features identified in the PK images. Figure [6] 
shows that the measured values depend on the type of so- 
lar feature, disk position, and filter bandpass. In particular, 
contrast values measured for each feature at all disk posi- 
tions in both PK-010 and PK-027 images are higher than 
those derived from PKR-011 observations. This fits in with 
the well-known observational evidence that the contrast of 
bright features is greatest at the line core. However, the 
contrast values measured in the PK-010 images, for the 
disk positions with [i < 0.5, overlap those obtained for the 
PK-027 observations, which sample more extensively the 
Ca II K line wings. 
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We found that the median of the measured contrast 
values, for a given feature over each \x range, increases to- 
wards the limb for all the features identified in the PK-027 
images (by 3% for the network, by up to 24% for the plage). 
The same quantities evaluated for all the features identified 
in PKR-011 observations appear to be constant across the 
disk, within the scatter of measurements. This is also true 
for the results obtained in the PK-010 images for plage and 
bright plage features (associated with reference models P 
and H, respectively). On the other hand, the network and 
enhanced network contrast measured in the same images 
(associated with reference models D and F, respectively) 
show a clear monotonic decrease toward the limb. This de- 
crease is likely due to both the limited spatial resolution 
of analysed images and a foreshortening near the limb. We 
find that the maximum value of contrasts measured for the 
same features and bandpass over each \x range is constant 
across the disk, within the scatter of the measurement un- 
certainty. Figure [7] shows the CLVs of maximum contrast 
values measured for the various disk features and band- 
passes. These results agrees qualitatively with published 
CLV measurement s of bright features from full-disk contin- 
uum observations (|Ermolli et al.ll200"7l ICriscuoli fc Ermollil 
2008, and references therein). 

Figures [6] and [7] also compare the measured and syn- 
thetic contrasts of models D, F, H, and P. Figure [B] (all 
panels) shows that the median of the contrast values mea- 
sured at the disk center are, on average, a factor » 1.7 lower 
than those derived from PRD using the corresponding at- 
mospheric models; the discrepancy increases for plage re- 
gions (associated with model H) identified in all the images 
and for all the features in PKR-011. Figure [5] (top and bot- 
tom panels) also indicates that the median contrast values 
deduced from PK-027 and PKR-011 images are in general 
poorer agreement with model predictions than the PK-010 
data. 

We also find that the values obtained with CRD for the 
models D to H, which are representative of network to plage 
regions, at disk center, on average, are only < 2% higher 
than those obtained with PRD on the same models, but 
the same quantity evaluated for model P, associated with 
bright plages, is « 29%, 20%, and 6% higher in PK-010, 
PK-027, and PKR-011 images, respectively. When moving 
toward the limb, the discrepancy between measured and 
modeled values slightly decreases for PRD computations, 
while it increases for CRD calculations. We find that the 
maximum contrast values measured for most of the solar 
features, filter bandpasses, and disk positions agree in gen- 
eral more closely with the outcomes of the spectral synthesis 
than the median contrast measurements. In particular, Fig. 
[7] (top and middle panels) shows that the CLV of the max- 
imum contrast values measured for plage and bright plage 
features identified in PK-010 images is very satisfactorily 
reproduced by the CLV of contrast values obtained with 
PRD for the corresponding reference models H and P, re- 
spectively. This is also true for the same features identified 
in PK-027 images, and for plages in PKR-011 observations, 
but with lower fidelity. 

5. Discussion 

The PSPT observations analysed in this study 
are characterized by high photometric precision 
(jRast. Ortiz and Meisne"rl 120081: iGoldbaum et al.l I2009D 



and a nearly simultaneous sampling of the solar atmo- 
sphere in various spectral bands. However, they have 
only moderate spatial resolution and miss the small- 
scale chromospher ic structure found by high-resolution 
observations (e.g. IWoger et al.l 20061 : Rezaei et"al~l 120081 : 
iCauzzi et al.l l200l iScharmer et al.l 12008ft . Thus. PSPT 
images are ideal for studying large-scale phenomena in 
the solar atmosphere, as in the present study, in which 
temporally and spatially averaged intensity measurements 
are compared with atmosphere models, but are less suited 
to study of the local dynamics, which tend to be averaged 
at these pixel resolutions. A potential weakness of our 
study is that PK profiles are up to about ten times 
wider than the bandpasses typically used for Ca II K 
obse rvations, e.g. of Lyot-type filters and spec trographs 
(e.g iKentischer et al.l 120081: lErmolli et all [2009;), making 
quantitative comparison with results from narrower band 
observations difficult. On the other hand, we note that 
the bandpass of the PK-010 images is similar to that 
of the Ca II K high-re solution observations from the 
Swedish Solar Telescop e (|Scharmc r et al.l [20031) analysed 
e.g. by iPietarila et all (|2009l ). while the bandwidth of 
PK-027 data is close to that of the Ca II H high-resolution 
imag es, from the SOT/BFI instrument (|Tsuneta et al.l 
| 2008f) onboard the Hinod e spacecraft, studied e.g. by 
lLawrence fc Cadavidl (|2010l ). 

To assess the importance of the filter width, we compare 
our measurements with results available in the literature 
for observations of the QS taken with spectrographs at the 
Ca II K line core, i.e. close to the spectral range sampled 
with the PK, but with a narrower bandwidth than that of 
PK. Figure [3] shows that our CLV measurements from PK- 
010 imag es are i n strik ing agreement with the results pre- 
sented bv lZirkerl (|1968t ). Moreover, our measurements from 
PK-027 and PKR-011 images, i.e. those obtained at the line 
core with the widest bandpass a nd in the red line w i ng, re - 
produce th e resul ts pu blished by White fc Suemotol ([1968), 
IWittmannl (|1976l ). and iGrigorveva et al.l (|2009h fairly well. 
Thus, the intensity CLV of QS appears to be relatively in- 
sensitive to the exact width or central wavelength of the 
Ca II K filter. However, it also implies that the mea sure- 
ments bvlWhite fc Suemotol (|l968l) . IWittmannl (|l976l ). and 
IGrigorveva et al. I (120091) were derived from data obtained 
in different observational conditions (e.g. spectral sam- 
pling, spatial resolution, seeing) and after different types of 
process ing (e.g . data calibration, straylight compensation) 
than in lZirkerl (| 19681 ). 

The close agreement found between our measurements 
and published results indicates that any systematic errors 
affecting our QS measurements, originating in the data 
analysed or caused by image processing applied, do not ex- 
ceed the scatter in previous measurements. We expect that 
this also holds for the measurements of the bright features 
identified at the disk center. However, at the limb, our mea- 
surements do not account for the dependence of feature size 
on both atmosphere height and disk position. 

To investigate the effects of feature misclassification on 
our measurements, we applied a second independent de- 
composi tion method to the data , specifically that devel- 
oped bv lNesme-Ribes et al.l ()1996f ). A direct comparison of 
the results obtained with the two methods is hampered by 
the different number of feature classes. However, the results 
obtained from this test support those presented in Sect. 4. 
The contrast values measured for the features identified by 
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the test-processing are, on average, lower than those pre- 
sented in Sect. 4 and in poorer agreement with the results 
derived from the synthesis at all disk positions. 

In the present study, we compared PSPT measure- 
ments with the outcome of spectral synthesis performed 
with a well-validated code on the FC09 atmosphere mod- 
els. We are aware that none of these models can real- 
istically reproduce the very inhomogeneous atmosphere 
depicted by high-resolution observations. Some spatial 
and temporal scales of atmosphere patterns displayed by 
high-resolution observations, as well as some amplitudes 
of temperature fluctuations, can be partly reproduced 
only by hydrodynamic models (e.g. ICarlsson fc Steinl 1995 , 



19971: IWedemever-Bohm et all 120041: iHansteen et alj |2007: 
Leenaarts et al. 20091 ). However, other temporally and spa- 
tially averaged properties of the solar atmo sphere deduced 
from o bservations, e.g. the ones discussed bv lFontenla et all 
(|2009fl . are reproduced fairly precisely by the static and 



one-dimensional model s, such as those con sidered in our 
study. As discussed by iRezaei et al.l (|2008l ). these models 
can be employed to derive a number of diagnostics (e.g. 
response functions to temperature perturbations) and to 
represent temperature and density variations as a function 
of height in the solar atmosphere, at least across the spec- 
tral range and at the temporal and spatial resolutions of 
observations analysed for our study. 

We find that the agreement between measured and mod- 
eled results for the intensity CLV of QS varies with the re- 
distribution formalism used in the spectral synthesis. In 
particular and somewhat surprisingly, the measurements 
are, on average, in slightly closer agreement with CRD com- 
putations than the PRD calculations. This may be because 
the atmosphere model s utilized in this stu dy were com- 
puted assuming CRD (jFontenla et al.l [20091 . If so, CRD, 
which results in enhanced emission towards the limb, com- 
pensates for the steep CLV that would result from the low 
density chromospheric temperature if the more physically 
realistic formulation of PRD were included. Therefore, the 
reference model needs to be modified to fit the data using 
the PRD formulation. 

To investigate the effects of adopting different reference 
models on the results obtained, we compared measured val- 
ues for the various disk features and bandpasses with the 
outcomes of the synthesis performed with PRD on earlier 
sets of semi-empirical models. In particular, we analysed 
the VAL3 and FAL3 series, and their subsequent upgrade 
FA06. Figure [9] shows the ratio of the CLV of median QS 
intensity values measured in PK images to the CLV de- 
rived from PRD computations for the quiet Sun models 
of the series indicated above. For all bandpasses, the in- 
tensity measurements are best reproduced by the synthe- 
sis performed on the QS model of FAL3 series (FAL3-C). 
For some bandpasses, synthesis results for the QS models 
of both VAL3 and FA06 series (VAL3-C and FA06-C, re- 
spectively) reproduce measured values more closely than 
achieved with the QS atmosphere (model B) of FC09. In 
addition, computations using the bright feature models of 
the FAL3 set (FAL3-F and FAL3-P) reproduce the contrast 
measurements for plage and bright plage features (associ- 
ated with models H and P of FC09, respectively) in the 
PK-010 images, within the standard deviation of measured 
values. In particular, Fig. [10] shows that computations us- 
ing the FAL3-F model reproduce the contrast measured for 
plage features in all the PK images, within the standard de- 



viation of measured values over many of the disk positions 
analysed. Figure [10] also shows that the same measurements 
are reproduced by the computations using model F of FC09 
more closely than those derived from model H of the same 
set, which implies that the correspondence between models 
and associated features could be improved yet further. 

These findings both confi r m and extend the results pre- 
sented by iGrigorveva et al.l ([2009). These authors com- 
pared measurements from spectroscopic observations of 
QS taken in the spectral range containing the Ca II H 
and K lines, with computations performed with RH and 
PRD using various atmosphere models. In particular, 
they analysed the VAL3-C and F AL3-C models, as we ll 
as the QS model pre sented by iFontenla et al.l (|2007l ). 
IGrigorveva et al.l (|2009h found that this latter model, which 
is quite similar to model B of the FC09 set utilized here, 
most closely reproduces their intensity CLV measurements 
in the Ca II K line wings, while at the line core the syn- 
thesis results for this model are in poorer agreement with 
measurements than results for earlier models. We note, 
however, that ou r stud y differs from that carried out by 
IGrigorveva et al.l (|2009h in several aspects. In particular, 
we analysed observations of various solar features, mostly 
observed in the Ca II K line core, while the other authors 
restricted their study to QS mainly observed in the line 
wings. Moreover, we analysed the most recent set of atmo- 
sphere models presented in the literature, in addition to 
the earlier series, while only the latter were considered by 
IGrigorveva et all (j2009ft . Finally, we utilized both the PRD 
and CRD approximations of the spectral synthesis. This al- 
lowed us to investigate the effects of code approximations 
and adopted models. 

Discussing the results obtained, we also consider that 
the PK images are contaminated by low levels of stray- 
light. By checking the falloff of intensity just outside the 
solar limb, we found that the images analysed in this study 
show the same contamina tion level as that estimated by 
Criscuoli k Ermolli (2008) in Ca II K observations from 
the Rome-PSPT. In particular, we found that the ratio 
of azimuthally averaged intensity values measured at 1.2 
solar-disk radii to the intensity measured at disk center is 
0.0167±0.0005, 0.0150±0.0005, and 0.0120±0.0005 for PK- 
01 PK-027, and PKR-011 re spectively. Therefore, follow- 
ing ICriscuoli &; Ermollil (|2008| ) we expect an increase of the 
order of 30% in feature contrasts, if they are measured in 
PK images restored for straylight degradation. 

To evaluate the effects of straylight degradation on the 
comparisons presented above, we developed two test anal- 
yses. Firstly, we studied straylight degradation effects on 
synthetic solar images. In particular, we applied a model 
point spread function (PSF) to synthetic images derived 
from our calculations with RH and the FC09 models. We 
measured the decrease in the intensity and contrast values 
derived from these synthetic images, hereafter synthetic at- 
mospheres. This provided a rough indication of the influ- 
ence of straylight contamination on measured contrast val- 
ues, by avoiding photometric uncertainties introduced by 
the iterative process of image-restoring methods. We also 
repeated our measurements for all the features considered 
in this study in PK images restored from straylight degra- 
dation. 

The details of the test analyses we performed are given 
in the appendix. To summarize, we found that the convolu- 
tion of the model PSF with homogeneous atmospheres, i.e., 
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synthetic atmospheres composed of a single feature at all 
disk positions, only leads to a significant decrease in inten- 
sity and contrast values modeled for the brightest features 
close to the solar limb. On the other hand, convolution of 
inhomogeneous atmospheres with the same PSF, i.e., syn- 
thetic atmospheres composed of quiet Sun and a number of 
bright features at various fj, positions, leads to a significant 
decrease in contrast values estimated for all the features 
and disk positions. In particular, convolution of the PSF 
with synthetic atmospheres composed of 10 small features, 
each covering 0.02 solar-disk radii, on average, leads to a 
decrease by nearly a factor two in the modeled contrast val- 
ues. This value is similar to the mean discrepancy between 
measured and modeled results for the various bright fea- 
tures presented in Sect. 4. In agreement, image restoration 
from straylight degradation leads to an average increase 
of « 18% over all the solar features and a slightly closer 
agreement between measured and modeled values. Model 
predictions and contrast values measured in restored im- 
ages at disk positions /i > 0.2 differ by a factor ranging 
from 0.9 to 2.7 (mean is « 1.6), depending on the solar fea- 
ture and filter bandpass, this compared to a factor of from 
1.1 to 3.1 (mean of 1.9) for un- restored images. 

While the results of our test analyses indicate that mod- 
eled and measured values are in slightly closer agreement 
after taking into account the effects of straylight contamina- 
tion, the model predictions are still approximately a factor 
1.45 higher than measurements obtained for all the solar 
features at disk center in restored images. 



6. Conclusions 

We have presented a study of the radiative emission of solar 
features across the Ca II K spectral range. We have anal- 
ysed moderate-resolution PSPT observations taken with in- 
terference filters that sample the Ca II K range with dif- 
ferent bandpasses. We have compared the results of CLV 
measurements for various disk features with the outputs of 
spectral synthesis performed with a well- validated radiative 
code and the most recent set of semi-empirical atmosphere 
models presented in the literature, as well as with earlier 
similar models. This led to the following results: 

(1) In general agreement with results presented in the 
literature, we found that the CLV of intensity values mea- 
sured in the Ca II K range for quiet Sun regions depends 
on both the filter bandpass and image quality. However, 
the measured dependence is quite small. The quiet Sun's 
CLV obtained in this study from analysis of the PK im- 
ages is in striking agreement with published results from 
observations taken with spectrographs. 

(2) The intensity response function for temperature per- 
turbations derived for all the filters considered in this study 
depends only slightly on the filter bandpass and reference 
atmosphere; the most sensitive to higher atmospheric layers 
are the ones derived for the filters with the narrower-band 
and the reference atmosphere of a bright plage (model P). 
This model, however, differs from the other models adopted 
in this study, not only in terms of its higher chromospheric 
temperature, but also its smaller value of micro-turbulent 
velocity. We note that even for the narrowest filter, PK- 
010, 58%-82% of the contribution at disk center is from 
atmosphere heights below 500 km, with the fraction being 
lower for brighter features, while the same quantity evalu- 



ated for PK-027, which has a bandwidth similar to that of 
the SOT/BFI Ca II H filter onboard Hinode, is 84%-94%. 

(3) RH calculations with the CRD ap proximation of 
the q uiet Sun atmosphere (model B) of iFontenla et al.l 
(|2009f ). on average, reproduce measurement results of inten- 
sity CLV for quiet Sun regions slightly more closely than 
PRD computations. However, RH calculations with PRD 
for earlier models of the quiet Sun atmosphere reproduce 
measured quantities with a similar accuracy. 

(4) The median contrast values measured for most iden- 
tified features, bandpasses, and disk positions are, on av- 
erage, a factor rs 2 lower than the contrast values derived 
from RH for the reference models considered in this study. 
Computations in PRD display a closer agreement with the 
data, which, however, is still unsati sfactory. Calculations 
for older sets of bright feat ure models (jVernazza et alJl98lt 
IFontenla et al.ll99 3. 2006) fit measurements from plage and 
bright plage regions as well as those derived from the cor- 
responding new atmosphere models. 

(5) For most of the identified features, disk positions, 
and filter bandpasses analysed in this study, the maximum 
contrast values measured in PK images are in general closer 
agreement with model predictions than for median contrast 
measurements . 

(6) The mean discrepancy between modeled and mea- 
sured values reduces, on average, by about 12% after taking 
into account straylight degradation effects on PSPT obser- 
vations. 

In conclusion, the success of CRD computations of 
the new atmospheres in reproducing some of the quan- 
tities analysed in this study does not prove conclusively 
that the radiative code used in this approximation or the 
adopted models realistically reproduce the observed quan- 
tities. However, it supports the combined use of this for- 
malism and models, as carried out for the spectral mod- 
eling of disk featur es over wide wavelength ranges (e.g. 
IFontenla et al.l[2009] ). On the other hand, the synthesis per- 
formed with PRD, on average, produces radiative estimates 
that are in closer agreement with measurements achieved 
in this study from bright disk features than those obtained 
with the less accurate CRD computation, in particular for 
the observations taken at the line core. This applies equally- 
well to synthesis results derived using either the most recent 
or the earlier atmosphere models, thus supporting simula- 
tions performed with PRD on earlier atmosphere sets for 
diag nostics across narrow spectral ranges of the Ca II K line 
(e.g. iReardon. Uitenbroek fc Cauzzl 120091 : iPietarila et al.l 
2009). The results derived from our investigation also in- 
dicate that the neglect of PRD, to some extent, can be 
compensated for by a change in the reference models. This 
is very similar to the ambi guous situation reported by 
iRutten and Kostikl (fl982h and lShchukina fc Truiillo Buenol 
(|200lD in the modeling of Fe I lines in the photosphere, 
where LTE and NLTE calculations can equally well fit the 
line profiles, as long as the temperature gradient in the at- 
mosphere is adjusted separately in each case. 
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Appendix A: Straylight contamination 

To evaluate the effects of straylight degradation on the compar- 
isons between modeled and measured quantities presented in Sect. 
4, we developed two test analyses. Firstly, we measured the effects of 
the straylight contamination induced by a model point spread func- 
tion (PSF) on synthetic solar images. The PSF used in our anal- 
ysis follows Criscuoli & Ermolli (2008); it includes one Lorentzian 
and three Gaussian components, which are defined by seven param- 
eters. We fixed six param eters of the PSF components derived by 
Criscuoli & Ermolli (2008), while setting the remaining parameter to 
reproduce the average value of the straylight level measured in PK 
images at 1.2 solar-disk radii. This parameter describes the weight of 
the Lorentzian relative to the Gaussian components. 

Figurc [A,ll shows the model PSF utilized in our analysis. We con- 
volved synthetic images constructed from intensities calculated with 
the one-dimensional FC09 models for the various disk features and 
bandpasses analysed in our study with the same model PSF. We con- 
structed two sets of synthetic atmospheres. The first set is composed 
of homogeneous atmospheres. The second set contains inhomogeneous 
atmospheres with two components, each atmosphere displaying quiet 
Sun and a number of features of the same type at various fj, positions. 
We describe in the following the results obtained by assuming 10 fea- 
tures, each covering 0.02 solar-disk radii at all fi positions, i.e. fore- 
shortening was not taken into account. For the various FC09 models 
at the PK-010 bandpass, Figs. lA~2l and lA~3l show the intensity values 
computed at various fj, positions in the original synthetic atmospheres 
and the atmospheres convolved by the model PSF. 

We found that the convolution of the model PSF with synthetic 
homogeneous atmospheres leads to a significant decrease only in the 
intensity and contrast values modeled for the brightest disk features 
close to the solar limb. Figure P\..4l shows that the discrepancy between 
measured and modeled CLV for QS increases when moving toward the 
limb if straylight contamination of PK images is taken into account, 
except when comparing the CLV measured for PK-027 and the one 
derived from CRD computations. Moreover, we found that PK-027 
measurements are in close agreement with the CLV obtained from the 
inhomogeneous atmospheres containing network features after their 
straylight contamination. 

On the other hand, we found that the convolution of synthetic 
inhomogeneous atmospheres with the same PSF leads to a significant 
decrease in contrast values estimated for all the features and disk 
positions. This is illustrated by Figs . I A . 5l and IAT61 which compare the 
results presented in Sect. 4, and the CLV derived from the convolution 
of the synthetic atmospheres with the model PSF. The contrast values 
estimated in the synthetic atmospheres convolved by the PSF are by 
up to nearly a factor two lower than those obtained in the original 
RH computations. Therefore, these results suggest a closer general 
agreement between median contrast values measured in PK and the 
contrast values derived from the RH synthesis than presented in Sect. 
4, if straylight contamination of small-scale disk features is taken into 
account. 

Our second test analysis was performed on PK images restored for 
straylight contamination by the deconvolution of original images with 
a model PSF derived from t he analysed data. We applie d the method 
and algorithm developed by Criscuoli & Ermolli (2008) to the whole 
sample of PK good images analysed in Sect. 4. 

Figures lA.71 and I A. 81 compare model results presented in Sect. 4 
and the CLV of contrast values measured for the various solar features 
on images corrected for straylight degradation. We found that image 
restoration leads to a slight increase in the average contrast values 
measured for all solar features with respect to the values derived from 
un-restored images, in particular for the features identified in PK-027 
and PK-010 images. In addition, we found an increase in the maxi- 
mum contrast values measured for all the features, filter bandpasses, 
and disk positions, in particular for the network features in PK-027 
images. We also found slightly closer agreement between model pre- 
dictions and contrast values measured in restored images than derived 
from un-restored images. We found that model predictions and con- 
trast values measured at disk center (fj, > 0.9) in restored images 
indeed differ by a factor ranging from 0.81 to 2.17, depending on the 
solar feature and filter bandpass, while the same quantity evaluated on 
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un-restored images ranges from 0.98 to 2.58. Modeled and measured 
values differ, on average, by a factor of 1.45 after taking into account 
straylight effects, compared to a factor of 1.71 for measurements in 
un-restored images at the same disk positions. 



CN 
O 
I 

CL 



O 

o 



o 
I 

Q_ 



1.0 


model P 


0.0 


'SWs******* 


* * 


- 


<s ^ 












0.5 












— 




0.0 


- - *W»W»W»I»W»W I <M 1 III-*-* 

, , , 1 , , , 1 , — 




— 
— t 


-41 _ to ' 


— h 






0.2 


model F 








-or" 






0.0 


— , — , — , — i — , — , — , — i — , — 




— t 




— V 






0.1 


model D 

:/>■*<■ a g. o 













1 f 

/ 


0.0 












-■§-- 





.0 0.8 0.6 0.4 0.2 0.0 



2.0 
0.0 

1.0 
0.0 

0.5 

0.0 

0.1 
0.0 



model 


P 








H 1- 


— i — i — i — i — i — i — i — i- 




model 


H 






HlgHHHiifniijj 


H 1- 


SfflffljjflHl O-ffl-ffl-J-m- ur " ifl " 
— , — , — I — , — , — , — hi- 


A— i — 1 

— i — i — i — i — i — i — 


model 


F 






p.imn-nminihuiiwiip-aBgiuumBrnrfflTra m"~BS W "8" — m— -- — — ■ 


model 










Bfflffli 


"■""■^ m tu 


D ID | 



1.0 0.8 0.6 0.4 0.2 0.0 





0.5 




0.0 


o 




1 


0.4 


cr 




(PK 


0.0 


o 


0.2 






c 


0.0 


o 




O 






0.1 




0.0 





















m^jgmy-^jw *^ ^ 36 * W m 








i S — »— 




model 


H 1 1 1 1 1 1— 

H 


— • 1 








1 y 










«-«-!— 








— i — i — i — i — i — i — i— 


— ' h 


-*- 




( 1 — 




model 


F 










/ 

_ / 








-K- 


i- 


~*- 

1 1 — 






— i — i — i — i — i — i — i— 


*-* 
— ' h 








model 


D 










/ ■ 








-*- 


* - * - 


— -a - 





1 .0 



0.6 0.4 



0.2 



0.0 



Fig. 6. CLV of median contrast values measured (symbols) 
for various disk features identified in PK-027 (top panel), 
PK-010 (middle panel), and PKR-011 (bottom panel) im- 
ages. For each bandpass, the various sub-panels show mea- 
surement results for solar features ordered by decreasing 
contrast, i.e. from the top panel results for bright plage, 
plage, enhanced network, and network regions, respectively. 
The error bars represent the standard deviation of measure- 
ments. For each disk feature and bandpass, the dot-dashed 
line indicates the average of values measured at disk posi- 
tions fj, >0.9, while dotted and dashed lines show the re- 
spective CLV derived from RH calculations with PRD and 
CRD and the reference model corresponding to the given 
feature, also indicated in the legend. In the top panel, the 
solid blue lines show the threshold values applied to the 
feature identification. 
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Fig. 7. CLV of maximum contrast values measured (sym- 
bols) for various disk features identified in PK-027 (top 
panel), PK-010 (middle panel), and PKR-011 (bottom 
panel) images. The error bars represent the standard de- 
viation of measurements. Legend as in Fig. [6] 
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Fig. 8. Comparison between intensity CLV of quiet Sun re- 
gions identified in PK images (lines, as described in the leg- 
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in the legend) from observations taken with spectrographs 
at the Ca II K line core. Details are given in Sect. 5. 
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Fig. 9. Ratio of the median CLV measured for quiet Sun 
regions in PK-027 (top panel), PK-010 (middle panel), 
and PKR-011 (bottom panel) images to the CLV derived 
from the synthesis performed with PRD using various at- 
mosphere models. Dotted, solid, dot-dashed, and dashed 
lines show results derived f rom t he quiet Sun mod e ls pre - 
sented bv IFontenla et all (120091) . IFontenla et all (|2006h . 
IVernazza et al. ~ (fl98lh . and Fontenla et al.l (jl993t) . respec- 
tively. 
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Fig. A. 2. Profiles of intensity values in the set of synthetic 
homogeneous atmospheres (solid lines) computed for the 
various disk features and PK-010 bandpass. In the various 
sub-panels, the different atmospheres are identified by col- 
ors and labels: yellow, green, blue, violet, and black show 
bright plage, plage, enhanced network, network, and quiet 
Sun atmospheres, associated with models P to B, respec- 
tively. Distance from disk center is given in disk radius 
units. The dashed lines show the profiles of intensity values 
in the various atmospheres after their convolution by the 
model PSF. 
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Fig. A. 3. Profiles of intensity values in the set of syn- 
thetic inhomogeneous atmospheres (solid lines) computed 
for the various disk features and PK-010 bandpass. Legend 
as given in Fig. IA.2I The dashed lines show the profiles of 
intensity values in the various atmospheres after their con- 
volution by the model PSF. Distance from disk center is 
given in disk radius units. 
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Fig. A. 4. Ratio of the median of quiet Sun limb darkening 
measured in PSPT images to the quiet Sun limb darken- 
ing derived from the synthesis performed for the PK-027 
(green, solid line), PK-010 (violet, dotted line), and PKR- 
01 1 (red, dashed line) filters by assuming either the PRD 
(top panel) or CRD (bottom panel) approximation. The 
lines indicate the results obtained by taking into account 
the CLV derived from the synthesis, while circles show the 
results derived by convolving the synthetic homogeneous 
quiet Sun atmosphere (associated with model B) by the 
model PSF. Asterisks show the ratio of the CLVs measured 
in PK to the ones derived from contaminated inhomoge- 
neous atmosphere composed of quiet Sun and network small 
features, associated with models B and D, respectively. 
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Fig. A. 5. CLV of median contrast values measured (sym- 
bols) for various disk features identified in PK-027 (top 
panel), PK-010 (middle panel), and PKR-01 1 (bottom 
panel) images un-restored for straylight degradation. For 
each bandpass, the various sub-panels show measurement 
results for solar features ordered by decreasing contrast, i.e. 
from the top panel results for bright plage, plage, enhanced 
network, and network regions, respectively. The error bars 
represent the standard deviation of measurements. For each 
disk feature and bandpass, the dot-dashed line indicates the 
average of values measured at disk positions fi >0.9, while 
dotted and dashed lines show the respective CLV derived 
from RH calculations with PRD and CRD and the reference 
model corresponding to the given feature, also indicated in 
the legend. Triangles and circles indicate the results ob- 
tained after convolution of the synthetic atmospheres by 
the model PSF, for homogeneous and inhomogeneous at- 
mospheres, respectively. In the top panel, the solid blue 
lines show the threshold values applied to the feature iden- 
tification. 13 
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Fig. A. 6. CLV of maximum contrast values measured 
(symbols) for various disk features (model P to model D) in 
PK-027 (top panel), PK-010 (middle panel), and PKR-011 
(bottom panel) images un-restored for straylight degrada- 
tion. Legend as given in Fig. IA.5I 



Fig. A. 7. CLV of median contrast values measured (sym- 
bols) for various disk features identified in PK-027 (top 
panel), PK-010 (middle panel), and PKR-011 (bottom 
panel) images restored from straylight degradation. For 
each bandpass, the various sub-panels show measurement 
results for solar features ordered by decreasing contrast, i.e. 
from the top panel results for bright plage, plage, enhanced 
network, and network regions, respectively. The error bars 
represent the standard deviation of measurements. For each 
disk feature and bandpass, the dot-dashed line indicates the 
average of values measured at disk positions fi >0.9, while 
dotted and dashed lines show the respective CLV derived 
from RH calculations with PRD and CRD and the reference 
model corresponding to the given feature, as also indicated 
in the legend. In the top panel, the solid blue lines show 
the threshold values applied to the feature identification. 
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Fig. A. 8. CLV of maximum contrast values measured 
(symbols) for various disk features (model P to model D) in 
PK-027 (top panel), PK-010 (middle panel), and PKR-011 
(bottom panel) images restored from straylight degrada- 
tion. For PK-027, measurement results for enhanced net- 
work and network features (associated with model F and 
D of Fontcnla et al. (2009), respectively) lie outside the 
plotted range, that is the range utilized for results from 
un-restored images. Legend as given in Fig. IA.7I 
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